Explosive percolation in scale-free networks 
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We study scale-free networks constructed via a cooperative Achlioptas growth process. Links 
between nodes are introduced in order to produce a scale-free graph with given exponent A for 
the degree distribution, but the choice of each new link depends on the mass of the clusters that 
this link will merge. Networks constructed via this biased procedure show a percolation transition 
which strongly differs from the one observed in standard percolation, where links are introduced just 
randomly. The different growth process leads to a phase transition with a non- vanishing percolation 
threshold already for A > A c ~ 2.2. More interestingly, the transition is continuous when A < 3 but 
becomes discontinuous when A > 3. This may have important consequences both for the structure 
of networks and for the dynamics of processes taking place on them. 

PACS numbers: 89.75.Hc, 05.45.Df 
Keywords: Networks, percolation 



The modern science of networks [l], 0, El has opened 
new perspectives in the study of complex systems. The 
simple graph representation, where the elementary units 
of a system become nodes and their mutual interactions 
links connecting the nodes pairwise, enables one to un- 
derstand a lot of properties about the structure and dy- 
namics of a system. In particular, the degree distribu- 
tion P(k), i. e. the probability distribution of the num- 
ber of neighbors k of a node, plays an important role. 
Real networks often display skewed degree distributions, 
where many nodes with low degree coexist with some 
nodes with high degree (hubs). The presence of the hubs 
is responsible for a number of striking properties, like 
a high resilience against random failures/attacks [4[ and 
the absence of an epidemic threshold [5[ . Resilience is de- 
termined by checking what is the fraction of nodes /links 
that need to be removed in order to split the network into 
a set of microscopic disjoint connected components. This 
is closely related to the process of percolation where 
one studies the conditions leading to the formation of a 
macroscopic (giant) component of the network. Here one 
starts from a set of nodes and no links; links are added 
randomly or according to a certain rule, until a giant 
component is formed. On networks having power law de- 
gree distributions (scale-free networks) with exponent A 
smaller than 3, the fraction of nodes/links to be removed 
from the graph for it to have no giant component tends 
to 1 in the limit of infinite network size [7( . In the spirit 
of percolation, and focusing on links, this can be read 
the other way around: a scale-free network with A < 3 is 
kept connected by a vanishing fraction of randomly cho- 
sen links, i.e. the percolation threshold is zero. For A > 3, 
instead, a finite threshold appears. Indeed, a giant com- 
ponent exists if the average number z 2 of next-to-nearest 
neighbors of a node exceeds the average number Z\ = (k) 
of its nearest neighbors Q . On networks without degree- 
degree correlations Z2 = (k 2 ) — (k), which diverges 
when the exponent A of P(k) is smaller than 3, whereas 
it is finite when A > 3. The divergence of the variance 



(k 2 ) is generally a sufficient condition to ensure the ab- 
sence of a percolation threshold on a scale-free graph, 
with or without degree-degree correlations [lfj , although 
in the case of large disassortativity a finite percolation 
threshold may emerge 11 1. 

In any case, whether there is a finite threshold or not, 
the percolation transition in networks is continuous: the 
order parameter, represented by the relative size of the gi- 
ant component with respect to the whole system, varies 
continuously from zero starting from the critical point. 
This is due to the fact that links are (usually) randomly 
placed on the network. Recent work by Achlioptas and 
coworkers has shown that, for networks similar to Erdos- 
Renyi random graphs [T3| , the percolation transition be- 
comes discontinuous (first-order) if links are placed ac- 
cording to special non-random rules [lj|. Such growth 
processes for graphs are meanwhile known as Achliop- 
tas processes, and the resulting connectedness transition 
as explosive percolation. Discontinuous transitions trig- 
gered by similar mechanisms were previously observed in 
the jam ming of information packets on communication 
networks 14| . In this letter we want to explore what hap- 
pens if one grows a scale-free network via an Achlioptas 
growth process. We will see that the resulting scenario 
is very different than in the case of ordinary percolation. 

Let us first define an Achlioptas growth process. The 
goal is to construct a random network of N nodes and 
given degree sequence {k\ , k%, . . . , kjy}. If links are placed 
randomly, the procedure can be carried out with the con- 
figuration model (l5j . Here instead, the criterion to add 
links is different. At the beginning of the algorithm 
(i.e., stage t = 0), we set fc s (0) = k s for each node 
s (the only condition needed is that ^ s k s should be 
an even number). The variables k s (t) act as a sort of 
counters: whenever a stub incident on node i is con- 
nected to another stub incident on node j, ki(t + 1) = 
ki(t) — 1 and kj(t + 1) = kj(t) — 1. The construc- 
tion proceeds until T — \ J2 S links have been drawn, 
which stands for k s (T) = , Vs (i.e. there are no 
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Figure 1: Scheme of the construction process of a network 
via an Achlioptas process with product rule (PR). Two pairs 
of stubs are taken at random (each pair is indicated by the 
dotted lines), and the products of the sizes of each pair of 
clusters merged by joining the stubs are computed. The stubs 
which are finally joined are those minimizing the product of 
the corresponding cluster sizes. In the case illustrated, one 
would join the nodes i\ and ji, which yield a smaller product 
cluster size than 12 and j'2 (2 ■ 5 = 10 versus 3 • 4 = 12). 



more stubs to be connected between node pairs). At 
each stage t of the growth, two pairs of vertices ji) 
and («2, 32) are selected as candidate links: these nodes 
are randomly selected among all vertices in the network 
with probabilities Pi t {t) = fc<i (<)/ E s p n {t) = 

M*)/E 8 M*)i PiM = M^/EsM*) andp j2 (f) = 
kj 2 (*)/ Ea k s (t) , respectively, which basically means that 
the candidate links (h,ji) and [i2,jz) are respectively 
selected with probabilities P(i 1 ,j 1 )(i) = PiA^Pjitt) and 
P(i2,j2)(t) = Pia(t)Ph(t)' I n or der to decide which of 
the two candidate links should be selected to become a 
real link to be added to the network, one computes the 
quantities L/^ j^t) = Mi t (t) Mj t (t) and Lu 2i j 2 \(t) — 
M i2 (t) Mj 2 (t) , expressing the product of the sizes of the 
clusters that the two selected links would merge (Fig.[T]). 
Finally, one draws the link for which the quantity L is 
lower. The former selection rule is called product rule 
(PR) . In principle other different reasonable criteria may 
be used instead of the PR: taking the sum instead of the 
product, maximizing instead of minimizing, etc.. Dur- 
ing the construction of the network, one should avoid 
the presence of multiple links (links connecting pairs 
of nodes already connected) and self-loops (links start- 
ing and ending at the same node). Scale-free networks 
may have a significant number of multiple links and self- 
loops 16], but in the transition regime we are interested 
in here they are essentially tree-like (most links have to 
be still placed), so multiple links and self- loops are very 
unlikely. In fact, we have verified that results do not 



change whether one allows or avoids them. 

On Erdos-Renyi graphs the process we have described 
generates a discontinuous percolation transition [l3j |. 
More recently, Ziff has studied the same process for bond 
percolation on two-dimensional square lattices [l7| , find- 
ing again a discontinuous transition. 

A natural parameter which allows to follow the con- 
struction of the network is p — t/T, which expresses the 
fraction of links added to the network during its growth. 
Following the construction of the network as a function 
of p allows to study the formation of the giant component 
and the associated percolation transition of the network. 
This technique allows to create the whole phase diagram 
of the transition through a single simulation [18j ]. 
Let us define as order parameter the percolation strength 
£(!) = A{W/N, where indicates the relative mass 

(i.e., number of nodes) belonging to the largest connected 
component in the network. If the transition is continuous 
(i.e., second-order), the theory of finite size scaling tells 
us that the percolation strength of a network composed 
of N nodes obeys the relation 



S {1) ^N-PI" f\{p- Pc ) N 1 ^ 



(1) 



where p c is the percolation threshold (in the limit of sys- 
tems of infinite size), j3 and v are critical exponents of 
the transition and F(-) is a universal function. Simi- 
lar laws of finite size scaling may be written for other 
observables. Here we consider the susceptibility x = 

N^J(SW 2 ) - (SW) 2 , which quantifies the amplitude of 
the fluctuations of the percolation strength. The suscep- 
tibility x obeys the relation 



X = N^G\(p-p c ) N 1 '" 



(2) 



where 7 is another critical exponent which characterizes 
the transition and G(-) is a universal function. The sus- 
ceptibility x is directly related to the order parameter 
S^l From the definition of x and the scaling behavior 
of at p c (Eq.[I]), we deduce that 7/1/ = 1 - (3/v. 
The susceptibility x ca n be used for the determination of 
the critical point p c . The percolation threshold p c (N) of 
a system of finite size N obeys the relation 



,(N) ^Pc + bN- 1 ^ 



(3) 



p c (N) can be determined by finding the value of p for 
which the absolute maximum of x occurs. Then a simple 
linear fit (based on the maximization of the Pearson's 
correlation coefficient) of p c {N) vs. N~ x / U allows to si- 
multaneously compute both values of p c and v. The same 
kind of analysis may be performed by determining p c {N) 
as the value of p at which one observes the absolute max- 
imum of (i.e., the relative size of the second largest 
component) [19(. In our numerical simulations, we find a 
perfect agreement between the two different approaches. 
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Figure 2: Explosive percolation transition in scale- free net- 
works. For A = 2.5 the transition is continuous. In (a) and 
(b) we show the percolation strengths corresponding to dif- 
ferent system sizes and their rescaling S^ 1 ' N 13 ^ , respectively. 
The validity of Eq.JTJ) can be proved by plotting S^N '" 
vs. (p — p c ) N 1/Il/ (c). The peak of the susceptibility x moves 
gradually towards p c as the system size increases (d). In- 
stead, for A = 3.5 the transition is discontinuous: percolation 
strengths corresponding to different system sizes do not have a 
scaling form (e). The location of the peaks of the susceptibil- 
ity is essentially the same for any system size (f). The network 
sizes go from 256000 to 16384000, via successive doublings. 



When the transition is discontinuous (i.e., first-order), fi- 
nite size scaling does not work. The scaling relations ([TJ) 
and @ trivially apply with (3/^ = and j/v = 1. The 
curves S 1 ' 1 ' vs. p corresponding to different system sizes 
do not scale and p c (N) approaches p c faster than as a 
power law in the limit of large N. 

We consider scale-free networks with degree exponent 
A [i.e., P(k) ~ fc~ A , where P(k) is the probability that 
a node has degree equal to k]. We examined two main 
scenarios, by setting the networks' cut-off (i.e., largest de- 
gree) equal to VN and to N 1 !^^ 1 ^. The results however 
do not qualitatively depend on this choice (the results 
shown refer to the cut-off V^V)- 

When a scale-free network is constructed via an 
Achlioptas growth process, the formation of the giant 
component is delayed. One needs to add a fraction of 
links much larger than in a standard random process be- 
fore seeing the emergence of the giant component. In- 
terestingly, for A < 3 it is already possible to measure 
a non-vanishing value of the percolation threshold. As 
an illustrative example, in Fig. [2] we show the behavior 
of the order parameter and the susceptibility x as a 
function of p. We see that for A = 2.5 the scenario is the 
one expected for a continuous transition, as confirmed 
by the scaling behavior of of Fig. ((2t), whereas for 
A = 3.5 the situation is different. 



Figure 3: Percolation threshold p c (N) as a function of the 
degree exponent A for various network sizes N. The black line 
represents the infinite size limit extrapolation of the critical 
threshold. This extrapolation is made by using Eq. <(3j for 
A < 3, while for A > 3 one uses the value of the plateau the 
pseudocritical points converge to. 



We have carried out a detailed finite size scaling analy- 
sis of the percolation transition in the range of exponents 
2 < A < 5. For each value of A we have determined the 
pseudocritical point at a given system size N and derived 
the infinite size limit of the threshold by using Eq. ([3]). 
In Fig. [3] we plot the lines of the pseudocritical points for 
various network sizes as a function of A. The black line 
indicates the extrapolation to the infinite size limit. The 
threshold is essentially zero up to A c ^ 2.2, and becomes 
non-zero for A > A c . From our analysis we cannot ex- 
clude that for 2 < A < A c the threshold is non-zero but 
very small; in order to clarify the situation one should 
use systems of orders of magnitude larger than the ones 
we studied, which lie already at the boundary of what 
one could do without using supercomputers. 

Interestingly, for A > 3 the pseudocritical point ap- 
proaches the actual threshold faster than as a power law, 
and the relation ((3]) does not hold, which hints to a first 
order phase transition. We have confirmed the result by 
performing the test suggested by Achlioptas et al. [HI]. 

In the region of A-values where we observe the second 
order phase transition we also computed the critical ex- 
ponents, by performing a finite size scaling analysis of the 
two main variables and \ & t the critical point, ac- 
cording to Eqs. dU) and (JSJ). We have used such analysis 
also to double-check in an independent way the extrap- 
olated values of the thresholds as a function of A, which 
we had previously obtained from the scaling of Eq. [3] 
the agreement is very good. The results are illustrated 
in Fig. |U We plot the values of the exponents' ratios 
"f/u and the sum (5/u + ^/u. We see that f3/v, "f/u 
are always in the range between and 1, but their val- 
ues depend on A. The sum (3/v + j/v is always 1 with 
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Figure 4: Critical exponents' ratios fl/v and y/v as a function 
of the degree exponent A. The insets show the scaling at p c 
of the quantities 5* (1) ~ N~ p/u (inset a) and X ~ N l/v (inset 
b) for A = 2.5, 3.0 and 3.5 (from bottom to top). 



good approximation, as expected. We also remark that 
around A c the exponents display a jump. This is due 
to the fact that the threshold goes to very small values 
for A < A c (consistent with zero), and finite size scaling 
cannot be accurate. For A > 3 the exponents take trivial 
values: (3 — 0, as the order parameter at criticality does 
not vanish in the infinite size limit; 7 = v, as the sus- 
ceptibility is an extensive variable, as it should be if the 
transition were discontinuous. The insets show the finite 
size scaling analysis on both exponents' ratios for three 
values of A. For we see that, while for A < 3 there is 
a clear power law scaling, as it should be for a continuous 
transition, for A > 3 there is a saturation. Similarly, for 
the susceptibility \, we see that the scaling is non-linear 
with N for A < 3, whereas for A > 3 it becomes linear, 
as it happens for extensive quantities. 

We have studied the percolation transition on static 
scale-free networks built with an Achlioptas process 
with product rule. We have found striking differences 
with standard percolation, from the existence of a finite 
threshold for A < 3 to the discontinuous character of the 
transition for A > 3. We stress that, since in an Achliop- 
tas process links are not placed completely at random, 
during the process the network generally has a different 
degree distribution, and only at the end of the process, 
when all links are pla ced, one restores the original im- 
posed distribution (20j. We have verified that the net- 
works at the percolation transition still have a power law 
degree distribution, but with a different exponent than 
the imposed one. In particular, we have verified that 
A = 2.2 corresponds to the effective exponent A' = 3. 
This may explain the existence of a finite threshold for 



A > 2.2 (it would correspond to A' > 3 for the actual 
networks at the threshold), but not the origin of the dis- 
continuous transition, which remains yet to be uncovered. 

Our findings show that the building mechanism of 
scale-free networks may strongly affect dynamic pro- 
cesses taking place on the network, along with structural 
features (e. g., resilience to failures/attacks), even if the 
degree distribution is predefined. So, very different phe- 
nomena can occur on networks with exactly the same 
degree distribution. The process we have studied here 
deserves further investigations, both from the numerical 
and the analytical point of view and it may reveal new ex- 
citing perspectives in the field of complex networks and in 
the theory of critical phenomena. Moreover, this finding 
may open new perspectives in other fields where networks 
are important, such as computer science and engineering. 
In particular, the issues of robustness and information 
transmission are inextricably linked to percolation. 

We are indebted to J. J. Ramasco for bringing this 
problem to our attention. S. F. gratefully acknowl- 
edges ICTeCollective, grant number 238597 of the Euro- 
pean Commission. At the moment of submission of this 
manuscript we have noticed a paper by Cho et al., just 
posted on the electronic archive |20j. The paper deals 
with the same problem, but the model used to build the 
network is not the same, which leads to significant dis- 
crepancies in the results. We apologize with Cho et al. 
for this unlucky and unwanted coincidence. 
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